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Recent HERA low Q^ data show that the logarithmic slope of the proton structure 
function ( q^^q2 ) is significantly different from perturbative QCD expectations for small 
values of Q^ at exeedingly small values of x. We show that shadowing ( screening ) 
corrections provide a natural explanation for this experimental observation. 
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Recent HERA data on Q^ and x dependence of the logarithimic Q^ derivative of 
the proton structure function F2(x,Q^), ^^^^2 , have become available [||]pl. This is 
shown in Fig.l where each point corresponds to a different values of Q^ and x. As seen, 
this observable rises steeply with decreasing x to values of x of approximately 10~^ with 
values of Q^ larger than a few GeV^ . However, in the exceedingly small x < 10~^ and low 
Q^ < 5 GeV^ values, the behaviour of the F2 slope is completely different - decreasing 
rapidly with x and/or Q^ getting smaller ( see Fig.l). This dramatic change of ^^^^2 
is not compatible with the prediction of perturbative QCD (pQCD ). This is shown in 
Fig.l where the HERA data ||l| is compared with the theoretical DGLAP expectations 
based on the GRV'94 parton distribution input [Q. Specifically, in the small x limit of 
the DGLAP equations we have 



(91n(QVA2) 977 



xG^^^^^(x,g^) , (1) 



where the expected rise of the F2 slope is associated with the logarithmic Q^ growth of 
xG{x, Q^) implied by the DGLAP equations in the small x limit. 

The new data pose a severe theoretical challenge as it requires a successful descrip- 
tion of DIS in the transition region characterized by intermediate distances. These dis- 
tances are smaller than the confinement scale ■^, where as(A^) ^ 1, but are still not 
small enough to justify a reliable pQCD calculation. For sufficiently small Q^ we expect 
non-perturbative contributions to be important in DIS. However, if the transition region 
corresponds to larger values of Q^ as it appears from Fig.l, we presume that these non- 
perturbative contributions to DIS will depend on several average characteristics such as 
vacuum expectation of the square of gluon tensor and the correlatrion length of two gluons 
within the hadron. These properties are reflected in the QCD Sum Rules approach |^. 

In order to assess the significance of the new data we examine both -^^hi and -F2. A 
change in the functional Q^ dependence of gf^ at small Q'^ values can be deduced from 
general arguments. Since at = ^"^02^ -^2 is finite ( non zero ) for real photoproduction 
( Q^ = and x =0 ), we conclude that at small Q^ values F2 oc Q^ ( see Ref. [^ for a 
more detail discussion). Accordingly, an eventual change in the functional x dependence 
of q^^q2 in the limit of exceedingly small x is expected. However, the particular x and 
Q^ location of this change, its dynamical mechanism and the relevance to pQCD are not 
apriori clear. We note that the F2 data in the HERA kinematic region at low x and 
Q^ > 1 GeV^ are well reproduced by the DGLAP evolution equations with input parton 
distributions 0| [0 HI starting from very low Q^ values. The importance of the new ^f^ 
data is that it opens a new window through which a transition region is observed in the 
Q^ range of 1 — 4 GeV"^ and x < 10~^. This transition is apparently not resolved through 



the study of F2 in the same x and Q^ domain. This transition is not predicted by the 
DGLAP evolution ( see Eq. ([l|) ) as can readily be seen in Fig.l. 

The objective of this letter is to show that shadowing ( screening ) corrections 
( SC ) to F2 account for the deviation from the DGLAP predictions and provide a natural 
explanation for the observed experimental phenomena. 
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Figure 1: The HERA data and GRV'94 prediction (triangles) on the F2{x, Q^) slope. The 
data are taken from Ref. 0]. 

There are two different types of SC that contribute to change of the F2 slope: 
(i) SC due to passage of the quark - antiquark pair through the nucleon, which lead to 



a more general equation for the F2 slope than Eq. ([1|) and (ii) SC to the gluon structure 
function in Eq. (|l]). We first discuss the quark- ant iquark sector. 

1. Closed formulae for the penetration of a gg - pair through the target in the Eikonal 
( Glauber ) approach were written many years ago |^ [|r^ and have been discussed in 
detail over the past few years [|ll[. For the F2 slope these formulae lead to 



dF^ix^Q^) Q' I ^^,^^ _ ^..(.,Q.;,.) ^ ^ ^2) 



where 



3Q2 



Here, we have used the main property of the DGLAP evolution equations, which allows 
us to factor out the impact parameter ( ht ) dependence from x and Q'^ ( see Ref. |]12[ ). 
Qq is the photon virtuality scale from which we apply pQCD for our calculation. T{ht) 
denotes the Fourier transform of the two gluon nonperturbative form factor of the target, 
which is independent of the incident particle. The relation between the profile T{ht) and 
the two gluon form factor reads 

nh) = - I e-^^'-^^^F{t)d\ (4) 

with t = —qf. Note, that factorization in Eq. (|^) is valid only for |t| < Qq. To simplify 
the calculations we approximate 

T{h) = -^e-i . (5) 

An impressive property of Eq. (|[) is the fact that the SC depend on the gluon structure 
function at short transverse distances r']_ = -^. This means that for the F2 slope we can 
calculate the SC with guaranteed theoretical accuracy in pQCD, while the SC for the 
deep inelastic structure function ( F2 ) which originate from large distances r]_ > -^ and 
are, thus, not well defined and could lead to errors in the calculations. 

The nonperturbative QCD ( npQCD ) information we need is given in Eq. (^ by 
the nonperturbative profile T{bt). In Fig. 2 we show the lowest order SC to F2 in Eq. (0) 
which are proportional to k^. In the additive quark model ( AQM ) we have two diagrams 
shown in Fig. 2 which are of order k"^. One can see that in AQM we have two scales for the 
integration over q±: the distance between two constituent quarks in a nucleon ( the first 
diagram in Fig.2 ) and the size of the constituent quark ( the second diagram in Fig. 2) 
( see Refs.||13|| ||l^ for details ). Eq. (H) is the simplest formula in which we can assume 



that two gluons inside a nucleon have no other correlation than their being confined in a 
nucleon with size R. A more general formula for the two radii model of the nucleon was 
obtained in Ref. H . 





Figure 2: The first order SC oc k^ for F2{x,Q'^). 
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xG^^^^^(x,g2)e ^ = Ki(x,Q2)e «! 
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(7) 



To evaluate the influence of the SC we introduce a damping factor ( D{k) ) 



ain(QVQd) 



^ j db^M^,Q^-bl) ' 



(8) 



where ^^^^2^,' :jj is calculated from Eq. (Q) for a one radius model of the nucleon and from 
Eq. (^ for a two radii model. The denominator is the first term of the F2 slope which 
corresponds to the DGLAP equations. Note, that in such a calculation for a two radii 
model R = Ri. 



Using this damping factor D{k,) we can write the F2 slope in the form 



dHQVQl 



DQif^) 



dHQVQl) 



(9) 



The calculated damping factor of Eq. (||) is plotted in Fig. 3 for a one radius model 
with R^ = 10GeV~^ |]I^ versus k ( upper curve ) and for two a radii model with two sets 
of radii: (i) Rj = lOGeV'^ and R^ = 3GeV-^ and (ii) Rj = 6GeV-^ and R^ = 2GeV-^ 
versus ki ( see Eq. (^ ). We note that two radii sets of curves are almost undistinguishable 
from one another as a function of Ki. 
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Figure 3: The damping factor Dq{k) versus k,{ki) = 0.03 + 0.05 n for a one radius model 
with R^ = 10GeV~'^ ( upper curve) and for a two radii model with two sets of radii: (i) 
Rl = lOGeV-^ and Rj = ^GeV-'^ and (ii) Rj = 6GeV-^ and R^ = 2GeV-\ 



One can see that the two radii model of the nucleon leads to sufficiently large SC 
which depend on the set of radii chosen. Note that the value of Ki is inversly proportional 
toRj. 

2. The Glauber ( Eikonal ) formula for the SC in the gluon structure function was 
obtained by Mueller [|l^ and discussed in details in Ref . [|^ . 

"^ dx' rQ' 

TT^ Jx X' 



xG"'^(x,Q^ 



-J ^ L dQ" dbf { 1 



-kg{x',Q'2;62) 



} 



(10) 



with KG{x',Q''^;bf) = ^K{x',Q''^;bf) where n is taken from Eq. (^). In the case of a 
two radii model the formula can be derived as a direct generalization of the procedure 



suggested in Ref. Wi 



xGix,Q') = -J'^ r dQ" Idhlil - e~-ori.<Q-:^\) + 

TX'^ Jx X' JQl J 



'IV. 



'^G2{^'^Q' jh) p-KGi(2^',Q'^;fe?) I 



t^Gl[x',Q'^■M) - ^G2{x',Q'^-bl) 



'^G2\X -iQ 'j^t, 



,KGi{x',Q'^;b^) - nG2{x',Q'';bl 



■g-«Gi(^',Q'';b?) _ g-«^G2(^',Q'2;b?) 



where kgi = 4^1 and kg2 = \i^2 where ni and ^2 are defined in Eq. (0). 
The gluon damping factor is defined as 

Ddx, Q') = . . .., J^y:^l^ , . ^. ... , (12) 
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/if^/S dQ'^JdblKGix',Q'^;bl) 
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where xG^^ is calculated from Eq. (|T0]) for a one radius model and from Eq. (|TT]) for a 
two radii model. For a two radii model kq = i^Gi in the dominator of Eq. (0). It is 
important to note that the Q'^ integration in Eq. (|ll|) and Eq. (|12]) spans all distances, 
including long distances dominated by nonperturbative dynamics. This is very different 
from Eq. (Q) and Eq. (j^) where the SC depend on the gluon density at short transverse 
distances. Since a theoretical approach to npQCD is still lacking, we eliminate the long 
distance contributions to Eq. ( jlQ]) and Eq. (|1T]) by imposing a low cutoff on the Q'^ 
integration. With this cutoff we neglect the contributions due to transverse distances 
r±_ > ^. This procedure makes our calculation of Dg{x^Q'^) less reliable than our 
calculation of Dq{x,Q^). To evaluate the errors due to this source we have calculated 
xG{x,Q'^) with two different cutoff values: Ql = OAGeV'^ and Ql = IGeV"^, using the 
GRV'94 parameterization |^ for the solution of the DGLAP evolution equations. The 
resulting gluon damping factors, Dg{x, Q"^) differ by about 10% which lends a reasonable 
credibility to our calculation. As we have noted in the introduction, we are of the opinion 
that the GRV'94 parameterization does not contain SC. This is due to the fact that most 
of the experimental sample used to fix the GRV'94 parameters has values of x > 10"'^, 
where we estimate the SC are very small. 

The formula which we use to compare with HERA experimental data is 

ff^..,«,..«^, «..,.,«=, f™2. ,13) 
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Figure 4-a: The F2 slope for SC in quark sector only. Triangles are the GRV'94 prediction. 
Stars are the result of SC calculations in the one radius model for a nucleon with R^ = 
lOGeV^, squares are for the SC in the two radii model with Rf = lOGeV^ and i?| = 3GeV^ 
and circles are for the SC in the two radii model with Rf = 6GeV^ and i?2 = 2GeV'^. 
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Figure 4-b: The F2 slope for SC both in quark and gluon sectors. Notations are the same 
as in Fig.4-a. 



where F^^^ is the deep inelastic structure function calculated in the DGLAP evolution 
approach with the GRV'94 parameterization. Our results compared with the experimental 
data |I| are shown in Fig.4 . Since we use the GRV'94 input our calculations can be carried 
out only for Q2 > QAGeV^. 

To summarize, the main points of this letter are: 

1) The deviation of q 1^(02/%) from the behaviour predicted by the DGLAP evolution 

in the small Q^ and exceedingly small x region is associated with SC applied both to the 
passage of quark - antiquark pair through the nucleon and the gluon density within the 
nucleon target. When SC are applied good agreement with the new HERA data, with 
Q2 > 0.8 GeV^ , is obtained. 

2) SC for the quark - antiquark interaction with the target nucleon is concentrated at 
short distances r_j_ ^ ^ and, therefore, can be calculated to a good degree of accuracy 
in pQCD. The effects of these SC are sufficiently large to account for most of the differ- 
ence between the DGLAP prediction and experimental data,, unlike the case of the deep 
inelastic structure function F2 p!3[] . 

3) The calculated SC for the gluon structure function are large and contain uncertain- 
ties due to long distance contributions which have not been included in the calculation. At 
present it is not possible to evaluate the errors which arise from nonperturbative contri- 
bution. We have checked the relative contribution coming from 3 GeV^ ^ r'j_ > 1 GeV^ 
by changing Qq - the lower limit of Q'^ integration in the calculation of xG^'"{x, Q^). The 
resulting change in Dg{x, Q"^) is not more than 10%. 

4) Our determination of the two radii of the nucleon rests on J/\I^ photo and DIS 
production data []T3| . The present analysis suggests that better data on the F2 logarithmic 



Q^ slope may provide an independent determination of these radii as well as additional 
knowledge on the role of long distance nonperturbative contributions to the SC. 
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